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CN Abstract 

Despite the no-go theorem [5] which rules out static stable solitons in Galileon 
theory, we propose a family of solitons that evade the theorem by traveling at the speed 
of light. These domain- wall-like solitons are stable under small fluctuations — analysis 
of perturbation shows neither ghost-like nor tachyon-like instabilities, and perturbative 
C ■ collision of these solitons suggests that they pass through each other asymptotically, 

which maybe an indication of the integrability of the theory itself. 

' 1 Introduction 

It is well-known that the special features of the DGP model [lj, which is a brane- 
world model aimed at modifying gravity in the infrared, are encoded in the dynamics of 
a scalar degree of freedom which has a Galilean symmetry — invariance under a shift of 
the field variable composed of a constant and a part linear in the spacetime coordinates i.e 
(j) — > (p+ao+a^x^.^. The Galilean model also arises in massive gravity theories where terms 
with Galilean invariance appear in the decoupling limit of massive graviton dynamics. [3] [I] 
It is an interesting question whether there are stable solitons in Galileon theory. In single 
field models, one could hope to find static solutions that stay in a local minima of the energy 
functional, and that , although they are not protected by topology, their classical stability 
is still guaranteed. However Endlich et al.[5] proved that in a generic single field Galileon 
model in any number of dimensions, there are no stable static solutions. Their proof rules 
out localized static lumps in 3+1 dimensional spacetime and, because the argument is valid 
also in 1+1 dimensional spacetime, it also rules out static domain walls. In this paper we 
show that in 1+1 dimensions there is a family of stable solutions that evade this theorem 
by traveling at the speed of light and then we study their collisions. There is a domain- 
wall counterpart for them in higher dimensions. The result shows that the Galileon model 
accommodates the existence of domain walls, but because they are moving at the speed of 
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light, they never look like static solutions. In section 2 we present the solutions and study 
their stability. In section 3 we find their energies and in section 4 we study the perturbative 
collision of these solitons and finally in section 5 we answer the possibility of existence of 
the solutions in higher dimensions. 



2 The Solutions and Their Stability 

The Galileon Lagrangian is: 

S£ = 8^4) + ad^d^&Dcj) + . . . . (1) 

Because the higher order Galileon terms vanish in 1+1 dimensions, we only present here the 
cubic interaction. The higher order terms can also be dropped for the domain-wall solitons 
in the 3+1 dimensional case. 
The equation of motion is: 

W + «[(n^) 2 - d^d^d^} = o . (2) 

Written in light-cone coordinate, with u = t + x and v = t — x, this equation simplifies 



to 



dudv 



f d 2 cj)\ d 2 (j)d 2 ^ 
\dudv J du 2 dv 2 



0. (3) 



It is obvious that if a function only depends on one of the light-cone coordinates, the 
equation is satisfied and we get the striking fact that any function of x ± t is a solution of 
([2]). Therefore we can take any localized lump as the initial condition and let it propagate 
at the speed of light. There is no rest frame in which the solution is static, therefore it does 
not contradict Derrick's theorem. 

Another way of seeing the existence of such solutions is to recall the fact that the Galileon 
theory has an interesting target-space diffeomorphism [6j[7j. The statement is that, for any 
solution (j)o(x, t) of the Galileon equation, an arbitrary smooth function / (4>o( x i t)) is also a 
solution of the Galileon equation if we only include the allowed highest order interaction in 
certain dimension[6j[7][9j.In two dimensions or for domain- walls in higher dimensions, the 
highest relevant non-linear term is the cubic one. Then one can see why these solutions 
exist: since 4>o(x,t) = x ± t satisfies the Galilean equation, therefore f(x ± t) is also a 
solution for an arbitrary function /. 

The main point of Derrick's theorem for the Galileon is ruling out the existence of the 
stable static solutions. Here we prove the perturbative stability of our non-static solutions 
against exponential growth of wave-like perturbations. We also eliminate the possibility 
of ghosts, waves with real frequency but the wrong sign in the action, by constraining the 
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overall sign of the action perturbation. This is the approach adopted by Nicolis et al.[2][8] 
in studying the stability of the DGP model. We apply the same approach here. 

Taking any classical solution (fio(x), perturbing around the solution and calculating the 
deviation of the action, we get 



5S 



d 4 xZ^ u du5(j)d u 



with 



Z" v = rf v + 2aU(j)Q^ u - 2ad^d\ 



(4) 



(5) 



Since 4>o is an extremum of the action, the term linear in 5(f) vanishes, leaving only terms 
quadratic and higher order in the perturbation. 

The stability is determined by the tensor Z^ v ' . If we look at Q as describing a free 
particle, Z^ v must have Minkowskian signature (+,—,—,—) to ensure that the particle 
has physical dynamics. Usually it is simpler to use the metric to turn Z^ v into Z%. The 
statement is then the Z% should be positive-definite. If all of the eigenvalues are positive 
after diagonalizing Zt, then the solution is stable. 

Written as a matrix, 



( 



\ 



1 - 2a<po -2a<p' \ 

2a</>' l + 2a<£g 

10 

1 



(6) 



First let's look at the right-moving solutions, i.e. 
In this case Zy simplifies to: 



f(x-t). 



Z» 



I 1 - 2a f" 2a f" \ 

-2a f" l + 2af" 

10 

1 



(7) 



The eigenvalue equation 



det(Z - XI) = (A - 1) 







(8) 



only admits 1 as an eigenvalue. The only condition we need for these solitons to be stable 
is the non-negativity of the eigenvalues, which is satisfied. This prevents the exponential 
growth of perturbations. However, there is a subtlety that we should address. The matrix 
Z^ is not diagonalizable and it can only be upper-triangularized, which can be a matter of 
concern. However, we can do the same analysis by perturbing the equation of motion and 
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study growth of the perturbations over time. 



( j ) (x,t) = f(x-t) + 5^(x,t) . (9) 

Plugging this into ^ gives: 

6(f>(x,t)(l-2af") -4af"5<j>'(x,t) - 5(f>"(x,t)(l + 2af") = . (10) 
Assuming a wave solution: 

5<f>(x,t) =Ae^ kx -^ . (11) 

leads to: 

2akf" ±k 

" = T^F' (12 > 

This u does not have an imaginary part and, based on the above arguments, we believe 
these solutions are locally stable. Here we treated /" as constant which means that the 
wavelength of perturbation is assumed to be small compared to the length scales of the 
variation of 4>(x). The same argument applies for the left moving solutions. For more 
details refer to M ■ 



3 Energy considerations 

The next step should be calculating the energy of these solitons. We can calculate 
the energy-momentum tensor of these objects, but before proceeding, we need to address 
subtleties regarding the higher order time derivative terms in the action. Expanding the 
Galileon Lagrangian in ([!]) gives 



dxdt 



<t> 2 - |V0| 2 + a U l <j> - <p 2 V z (j) - 0| V0r + V l <t>\V<t>\ 



(13) 



We assume that the solution is trivial at temporal and spatial infinities. The term (p 2 cp is a 
total time derivative and can be dropped. The <j> |Vi^| 2 term can be simplified 



dt 



|V0rJ -20(V0- V</> 
-V 



• V<A 



;,2 



2 V0 +0 2 v 2 < 



Dropping the total derivatives simplifies the action to 



(14) 



S 



dxdt 



\V(t)\ l - 2o^ 2 V 2 </> + aV^| V0f 



(15) 
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The conjugate momentum II and the Hamiltonian are 



n = S 4 = 2(1 - 2aV 2 

5<j) 



H 



dx 



L 



dx 



n 2 

4(1 - 2a0, a 



+ \<t>.x\ 2 - a6, xx \(f). x \ 2 



(16) 
(17) 



We use the right moving solution 4>(x,t) = f(x — t) to compute the Hamiltonian and get 



E[f] = J dxf'\2 - 3a/") 



(18) 



The second term is a total derivative and can be dropped and the energy simplifies to 



E[f] = 2 J dxf' 2 . 



(19) 



This energy is positive and remarkably does not depend on the interaction terms. The 
energy is exactly the same as if there were no higher derivative interactions in the Lagrangian 
in Q. The steeper the lump is, the larger the energy. The lump energy scales as the inverse 
square of its width. We can obtain the same conclusion by looking at the energy momentum 
tensor 



5S 



'- /it/ 



Si]' 



d^4>du(f> + an^d^dvcf) + ad\(f>d (f)d^d u <f> 



(20) 



For the right moving solution f(x — t), we have the stress tensor: 



- )1V 



f 2 -f' 2 
-r 2 r 2 



(21) 



Again we see that the energy momentum tensor is independent of the non-linear terms 
proportional to a and therefore satisfies all the energy conditions for free waves moving at 
speed of light. Despite the fact that in higher dimensions more interaction terms survive 
in Galileon lagrangian, they do not change this energy formula and it remains correct in 
higher dimensions. 



4 Perturbative collision of two solitons. 

One important aspect of studying solitons is their collision. If a right moving and a 
left moving soliton collide, the result could be highly nontrivial. Let's assume that initially 
there are two such lumps moving towards each other without any interference term, so that 

(t>{x,t) = f(x + t) + g(x-t) . (22) 
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This wave-function does not satisfy the equations of motion. The full solution should be of 
the form: 

(j)(x, t) = f(x + t) + g(x -t) + h(x, t) . (23) 

We try to solve this equation in a limit \h(x, t)\ <C \(f>(x, t)\ with boundary conditions 
limt^-oo \ h(x,t)\ = 0. The natural way to solve this equation is in the light-cone coordi- 
nates. Plugging (23) in ([2]), we obtain: 



h uv (l + 2ah uv ) = 2a(f" + h uu )(g" + h vv ) = . (24) 

If h and a are small, we can set up a perturbative expansion for this equation. Because 
a has dimensions of [Length] 3 , we need to be careful about the meaning of its smallness. 
The only relevant length-scale in this problem is the width of the two lumps. In this 
approximation, the smallness of a should be understood compared to the width of these 
lumps. If this criteria is not met, we cannot proceed much further. In this limit, we can 
iteratively integrate the equation of motion for h, assuming that it vanished far in the past. 
Expanding h : 

h = hP> + hW + ... . (25) 



The first order equation : 



whose solution is: 



h^\ uv = 2af"{u)g"{v) . (26) 



h<»(u, v) = 2af{u)g'(v) + 2af(u )g'(u ) - 2af{u )g'{v) - 2af(u)g'(v ) (27) 
+ h 1 {u ,v) + h 1 (u,v ) -h l {u Q ,v ) . (28) 

where uq and vq are the initial points for the integration. Choosing these points in the far 
past, i.e. uq,vo — > — oo, and assuming that / and g are localized lumps, it is easy to see 
that all the terms on the right except for the first one vanish. Using the first order equation 
for h , the second order equation is 

h&, uv = 2ah] uu g"{v) + 2ah] v J"(u) . (29) 

which gives 

h {2) (u, v) = 2a 2 f"{u)g ,2 {v) + 2a 2 f' 2 (u)g"(v) + Boundary terms . (30) 
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Using the fact that h vanished at far past, the boundary terms should vanish. Similarly the 
3rd order equation is: 



/i (3) , u „ = 2a[g"h 2 , uu +f"h 2 , vv +h 1 , uu h l , m -2ah l , uv f'g"] . (31) 

leading to: 

h®(u,v) = 4a 3 [f"'(u)g' 3 (v) + f\u)g"\v) + 6f' (u)f (u)g' (v)g" (v)] (32) 

/u pv 
f"(u) 2 du / g"(v) 2 dv + Boundary terms . 
-oo J — oo 

It is important to notice that if the width of the lumps is large compared to a, the higher- 
order terms in the perturbation, which contain higher derivatives, will be smaller . All the 
terms in the perturbation are interference terms between the right and left moving waves. 
Therefore, h is negligible except for during the time that the two waves pass each other. 
When they get far enough apart, the interference term goes back to zero. Therefore in the 
limit of small a there is no dissipation and the two solitons retain their shapes. In Fig|T] 
we showed the collision obtained by this perturbative method before, during, and after the 
collision of the two wave-packets of the form (a = 1) 

u 2 

f(u) = A[u(u - cj)i)(u - 4> 2 )(u - <f) 3 )(u - 4>±){u - cj) 5 ) + cf> 6 ] exp [ — — -A (33) 

width 

v 2 

g(v) = A[v{v + ct>i){v + fa)(v + <fo){y + fa)(v + S ) - 06 exp[ —J , (34) 

width 

where 

A = 1 width = 0.3 (f>i = -0.3 4> 2 = -0.1 
fa = 0.2 <p A = 0.3 <fc = 0.35 6 = 0.0015 . 

These graphs show that the two waves pass through each other perfectly and their shape 
and profile remains the same as before the collision after a long time. This can be an 
indication of the integrability of the Galilean Lagrangian as hinted by [6j. 



5 Conclusion 

With the above analysis, we have found a family of classical solitons in Galileon theory 
in (1+1) D. These solutions are locally stable and they have corresponding domain-wall-like 
structures in (3+1) dimensional spacetime. During collisions, these domain-walls would not 

1 When this paper was in preparation, we were notified by K.Hinterbichler that the same solutions and 
their stability were independently discovered by J.Evslin[§j, who studied closed time-like curves in this 
model. 
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Figure 1: Collision of the two solitons. The progression is from top left. As seen in these 
graphs, the two solitons retain their shape after they pass each other. 
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change shape or dissipate, at least in the case that the width of the walls are large compared 
to the coupling constant a. While there is a Derrick's theorem regarding the existence of 
stable static solutions in Galileon theory, we notice that stable "moving solutions" are not 
necessarily ruled out and this work provides an example. Regarding the Derrick's theorem, 
another kind of solutions that could not be ruled out are the gauged ones. Just as in the 
case of interactions via a potential, gauging the theory |12| may violate the Derrick-type 
arguments and allow the existence of stable static solutions. We also tried to find solutions 
analogue to a bubble moving at the speed of light in higher dimensions and could show the 
solutions are non-local and therefore not interesting. 
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